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Abstract. We describe theoretically the main characteristics of the steady state regime of a type II Optical 
Parametric Oscillator (OPO) containing a birefringent plate. In such a device the signal and idler waves 
are at the same time linearly coupled by the plate and nonlinearly coupled by the \^ crystal. This mixed 
coupling allows, in some well-defined range of the control parameters, a frequency degenerate operation 
as well as phase locking between the signal and idler modes. We describe here a complete model taking 
into account all possible effects in the system, i.e. arbitrary rotation of the waveplate, non perfect phase 
matching, ring and linear cavities. This model is able to explain the detailed features of the experiments 
performed with this system. 



PACS. 42.65.-k Nonlinear optics - 42.65.Yj Optical parametric oscillators and amplifiers - 42.60.Fc Mod- 
ulation, tuning, and mode locking - 42.25.Lc Birefringence 



1 Introduction 

In a type II OPO, signal and idler fields of crossed polar- 
izations are generated when the pump exceeds a certain 
threshold. Energy conservation requires that luq = c<Ji+a;2, 
where loq, oj\ and L02 are respectively the pump, signal 
and idler frequencies. The precise values of the signal and 
idler frequencies are set by the conditions of equal cav- 
ity detunings and minimum oscillation threshold, which 
depend on the value of the phase matching between the 
three waves, and on the vicinity of cavity resonances for 
the signal and idler modes. Theses frequencies are deter- 
mined unambiguously when one knows the values of two 
control parameters of the OPO, namely the crystal tem- 
perature (which sets the value of the different indices) and 
the cavity length (which determines the cavity resonance 
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conditions). Frequency degeneracy, i.e. io\ = L02 = ^o/2, 
occurs only accidentally since it corresponds to a single 
point in the parameter space. It cannot be achieved for 
a long time in real experimental conditions, as these pa- 
rameters drift in time. Furthermore, even when the device 
is actively stabilized on the frequency degeneracy work- 
ing point, the signal and idler fields still undergo a phase 
diffusion phenomenon, similar to the Schawlow-Townes ef- 
fect in a laser but acting on the difference between 
the phases of the signal and idler modes in the case of the 
type II OPO. As a result, the output field polarization 
direction slowly drifts with time. 

In the context of quantum information and generation 
of EPR correlated bright beams, where both amplitude 
and phase correlations are involved, phase locking is inter- 
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esting since it allows a much simpler measurement of am- 
plitude and phase quantum correlations between the sig- 
nal and idler beams |3j even above threshold : the measure- 
ment of intensity quantum correlations between the sig- 
nal and idler modes can be done even with non-frequency 
degenerate beams ^ but the measurement of phase cor- 
relations makes it necessary to use a local oscillator. Thus 
a phase reference is defined and signal and idler must be 
stable compared to this reference, which is not the case in 
a regular, above threshold, type II OPO. 

A few years ago, Wong et al. had the idea of achiev- 
ing the frequency degenerate operation at the output of a 
type-II OPO by introducing a linear coupling between the 
signal and idler fields. This coupling was made by way of 
a birefringent quarter-wave plate placed inside the OPO 
linear cavity which couples the two orthogonally polarized 
signal and idler waves. In this way, they generated intense 
and stable frequency degenerate signal and idler beams |o]- 
The theoretical model described in reference [B] was able 
to account for the main features of this phenomenon, but, 
for the sake of simplicity, it was made for a ring cavity, 
for a small angle between the crystal neutral axes and the 
birefringent plate neutral axes, and without any phase- 
shifts introduced by the reflection on the cavity mirrors 
or by non perfect phase matching. 

Most experiments use linear cavities, whereas most 
theoretical treatments assume ring cavities. For scalar fields 
there is almost no difference between the two configura- 
tions (the crystal in the ring cavity being taken twice as 
long as the linear cavity), if one neglects the cavity mirror 
differential phase-shifts. This is no longer the case when 
polarization effects are taken into account : in this case, 
a matrix formalism is needed, and the exact succession of 
the different elements in the cavity is now important, as 
they are described by non-commuting matrices. It is also 
interesting to examine the regime when the birefringent 
plate angle is not limited to small values. It seems also 
important to take into account the mirror phase shifts, 
which are known to induce a significant change in the 
phase matching between the three waves and consequently 



in the oscillation threshold of the linear cavity OPO 0. 
The purpose of the present paper is to introduce all these 
refinements in the theoretical model introduced in 6 , and 
also to discuss the properties of the phase-locked OPO 
in terms of the actual control parameters of the device, 
which are the cavity length and the crystal temperature. 
This paper is followed by a second one [HJ in which the 
quantum fluctuations and correlations between the signal 
and idler fields are determined and studied in the same 
configuration. 

In section Q, we introduce and describe the behavior 
of the different elements placed inside the OPO cavity. We 
then determine and discuss in section Q the steady-state 
regime in the ring cavity case. Finally, in section Q, we 
examine and discuss the steady state regime in the linear 
cavity case. 

2 Linear and nonlinear elements in the OPO 
cavity 

We consider here a x^ 2 -* crystal with a type II phase match- 
ing, which means that the signal and idler fields have 
orthogonal polarizations. The crystal length is / and its 
indices of refraction are ri\ and n 2 respectively for the sig- 
nal (ordinary) and idler (extraordinary) waves which are 
supposed to be frequency degenerate. The non degenerate 
case will be studied elsewhere [S]. 

Assuming a small variation of the various field ampli- 
tudes inside the nonlinear medium, which is quite reason- 
able in a c.w. OPO, one can solve in an approximate way 
the propagation equations inside the crystal, and obtain 
to the second order in the non-linearity, g : 

A (l) = MO) ~ ffexp (-i4M) sine (4^) ^(0)A 2 (0) 
_Z r (^i)(| Al (0)| 2 + |A 2 (0)| 2 )A o (0) 
A 1 (l) = A 1 (0)+ 5 exp(»^)sinc(^)A (0)A*(0) 

+ 4/(^)(|A (0)| 2 -|A 2 (0)| 2 )A 1 (0) 
A 2 (l) = A 2 (0) + geiq>(i^) sine (^)A (0) A* (0) 

+ ^/(^)(|A (0)| 2 ~|A 1 (0)| 2 )A 2 (0) 

(1) 
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Fig. 1. p is the angle between the crystal's axes ((Ci, C2), 
black lines) and the waveplate axes ((Wi, W2), grey lines). 

in which Aq is the envelope amplitude and the Ai, i = 
1, 2 are the envelope amplitudes of the interacting fields, 
assumed to be plane waves, along the crystal axes (C\ : 
ordinary wave, Ci : extraordinary wave); the envelopes are 
normalized in such a way that \Ai\ 2 gives the photon flow 
(photon. m _2 .s _1 ); g is the nonlinear coupling coefficient 
given by 

(2) 



(2) I hlU U!iUJ2 

V 2c i e n nin 2 



and f{x) = exp ^ x ^ (exp(ia:) — sinc(x)). The crystal input- 
output equations are then, when one expresses the pump 
field at the center of the crystal : 



A 1 (l) = A 1 (0) + g'A (i)Ai(0) 
A 2 (l) = A 2 (0) + g'A (l)At(0) 



(3) 



where g' = gexp sine The advantage of this 

expression is that it is valid to the second order in the 
non-linearity g' . 

The second element in the cavity is the birefringent 
wave plate. It has a thickness e and its indices of refraction 
are n e and nf at frequency luq/2 respectively for the slow 
and fast axes which make an angle p with the crystal axes 
(see fig. P). 

Its effect will be described in the Jones matrices for- 
malism ;9 in the nonlinear crystal axes basis. The trans- 
mission through the wave plate can be written as the ma- 
trix : 

a e 



M 



where 



n s + rif 



(4) 



(5) 
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Fig. 2. Set-up of the ring cavity; we consider only one direction 
of propagation. 

represents the mean index of refraction of the waveplate, 
and 



cos 



Zi0 

Ac/) 



i cos(2p) sin 



sin(2p) 



(f 



(6) 
(7) 



where A<p = k(n s — nf)e is the waveplate birefringent 
phase-shift. Let us set a = aoe^ where (cto,ip) £ R 2 . For 
the sake of simplicity, we will assume that this plate has 
no effect on the pump field polarization, i.e. acts as a A 
waveplate at the pump frequency. 

3 Ring cavity type II OPO 

We assume in this section that the cavity has a ring shape 
(fig. [21, and that the coupling mirror has large reflec- 
tion coefficients for signal and idler modes (ri and r 2 ). 
The moduli of the amplitude reflection coefficients of the 
coupling mirror are taken equal for the signal and idler 
modes : |n| = \r 2 \ — r — 1 — ft, with k <C 1, so that 
the transmission is \t\ ~ \2k. p. is the round-trip loss 
coefficient for the signal and idler waves (due to crystal 
absorption, surface scattering, other mirror finite trans- 
mission ...), assumed to be small. We define a generalized 
reflection coefficient : r' = r(l — /j,) w 1 — (/i + ft). We will 
call Ci and £2 the phase-shifts introduced by the reflection 
on the cavity mirrors for the signal and idler waves so that 
rj = r exp(i£j), j — 1, 2. In all the article, we do not take 
into account the resonance of the pump mode : all equa- 
tions are given for the pump field inside the crystal and 
we only calculate operating thresholds (not signal or idler 
intensities) normalized to the intracavity pump threshold 
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of the OPO without the waveplate (standard OPO thresh- 
old), (To. The free propagation length inside the cavity is 
denoted L. 

As the signal and idler fields are assumed to have the 
same frequency, the birefringent plate and the non-linear 
crystal couple the same fields, which are the ordinary and 
extraordinary waves at frequency ^ , and only three com- 
plex equations are needed to describe the system. From 
equations © and one readily obtains the following 
steady state equations for the field amplitudes A\ — j4i(0) 
and A 2 = A 2 (0) : 

Ax = r'aoe 1 ^- / 2 ^ {A x + g'A A* 2 ) + 
r'ee^+ e / 2 HA 2 +g'A Al) 

A 2 = r'a e^ +e / 2 -^ (A 2 + g' A A\) + 
r > ee Ks-e/2) ( v 4 1 + g'A A* 2 ) 

where 8 = f§ (f e + + L)+&±& is the mean round- 

trip phase-shift, and 9 = — n 2 )l + Ci — C2 is the bire- 

fringent phase-shift introduced by the non-linear crystal 
and by the mirrors. 

One immediately observes that these equations are not 
invariant under the gauge transformation A\ — > Aie lip , 
A 2 — > A 2 e~ lip , as is the case for the usual equations 
of a non-degenerate OPO without birefringent mixing. 
This implies that, unlike in the usual OPO, the phases 
of the signal and idler amplitudes solutions of equations 
©, when they exist, are perfectly determined : phase- 
locking has occurred between the two oscillating modes, 
and there is no phase diffusion effect. This phase- locking 
phenomenon is common to all linearly coupled oscillators 

EDI 

Since the effect of the different elements on the polar- 
ization is described by matrices which do not commute, 
one expects that the system depends on the plate position. 
However, it is straightforward to show that exchanging the 
positions of the waveplate and of the crystal amounts to 
a rotation of ir/2 of the crystal which is equivalent to ex- 
changing indices 1 and 2. This does not change the physics 
of the system so that we will place ourselves in the case 



where the waveplate is located after the crystal with re- 
spect to the input beam. 

Equations JSJl have been solved analytically in the small 
angle regime p <C 1 and for small cavity detunings and 
losses in ref 

We will present here the properties of the more com- 
plex analytical solutions obtained without any approxi- 
mations : we will not give the complicated expressions of 
the solutions, but instead give plots of the most striking 
results. The exact expressions for the different parameters 
in the case of a small angle are given in the appendix. 

The real and imaginary parts of the first two equations 
of (JSJ form a set of two linear equations for the amplitude 
and phase of the field enveloppes Ai, A 2 . Thus, one ob- 
tains a set of four linear equations with four variables. 
A non-zero solution of this systems exists only when the 
corresponding 4x4 determinant is zero. This condition 
gives a real equation for the system parameters, which 
is fulfilled only in a specific operating range, or locking 
zone, for the self-phase-locked OPO. In the locking zone, 
this equation has two real solutions for the intracavity 
pump intensity, corresponding to two possible regimes of 
the system In this paper, we will focus our attention 
to the regime of lower threshold. These solutions give the 
oscillation threshold for the intracavity pump power as a 
function of the crystal temperature, the cavity length and 
the waveplate angle. We define a as the ratio of the in- 
tracavity pump power to Co & n d o~ th as the ratio of the 
intracavity pump power threshold to o-q. If, for a given 
set of parameters, a is larger than a , one obtains fre- 
quency degenerate oscillation. One can thus plot the val- 
ues of cavity length and crystal temperature for which a th 
is smaller than a so that there is degenerate oscillation. 
Fig. (JSJ-a) displays the locking zones for two values of the 
wave plate angle p as a function of ST = T — T^ eg and 
5L = L — Ldeg, where Td eg is the temperature for which 
the exact frequency-degenerate operation occurs without 
any birefringent coupling and and Ld eg the correspond- 
ing cavity resonance length. The locking zone consists of 
two surfaces which overlap for small values of p. Fig.|(21b) 
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shows the cross section AA' of the locking zone for a given 
value of ST, that is a th as a function of SL. All curves in 
this paper are plotted in the case of KTP for which the in- 
dex of refraction vary with the following dependance jllj : 



dni 
~dT 



1.3xl0~ 5 K' 1 and 



dn 2 
~dT 



l.exKT 5 ^- 1 (9) 



5L (/jm) 
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Fig. 3. (a) : Locking zone as a function of the cavity length 
(SL) and of the crystal temperature (ST) for waveplate angle 
p = 1° (light grey) and p = 5° (dark grey), a = 3. (b) : a th as 
a function of the cavity length (SL) for ST = 0.2 K. For a — 3, 
this corresponds to the cross section AA' of the locking zone. 

A(f> = 7T. 

Fig. (JSJ-a) shows that the locking zone extension in- 
creases as a function of p. However, the minimum thresh- 
old does not increase with p and a threshold equal to the 
standard OPO threshold can always be found for ST = 0. 
For ST 0, the minimum threshold as a function of SL is 
no longer equal to one (see fig. El bottom). 



For a given value of p, the coupling parameter e = 
i sin( -4r) sin(2yo) is maximized for Acf) = ir that is for a A/2 
waveplate. As shown on fig. (@J, a different value for A<j> 
will reduce the locking zone extension but does not change 
the general shape. In order to maximize e and thus the 
locking zone extension, one can set A<f> — n and p = 45°. 
In this case, the locking zone is infinite, in practice only 
limited by the phase matching. 
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Fig. 4. Locking zone as a function of cavity length and crys- 
tal temperature. The thin dark grey line corresponds to a A/2 
waveplate while the light grey zone corresponds to a A/4 wave- 
plate. p — 5°, a = 2. 

As the locking zone depends on the temperature, it 
may be important to take into account the phase match- 
ing. However for small values of the waveplate angle, the 
locking zone extension in ST is small so that the effect of 
deltak ^ remains small. As p is increased, this effect 
becomes noticeable and limits effectively the extension of 
the locking zone to a zone ST w 10 K. We have plotted on 
fig. (J5J cr res , the threshold on resonance : it corresponds 
to the minimum value of a as a function of SL for a fixed 
value of ST. One notices on this figure that <r res is pe- 
riodic in ST if one does not take into account the phase 
matching : this is due to the periodicity in temperature of 
the crystal birefringence. When one takes into account the 
phase matching, this periodicity disappears (grey curve). 
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Fig. 5. Normalized threshold on resonance <r res as a function 
of the temperature for p — 30°. The black curve corresponds 
to the result obtained without taking into account the phase 
matching multiplied by 50 for readability of the figure. The 
grey curve is plotted taking into account the phase matching. 



crystal 



M, 



M, 



waveplate ' 



Signal/idler 



Fig. 6. Set-up of the linear cavity type II OPO 

For realistic parameters, such as R = 90% and a = 2, 
the transverse width of the locking zone is AL ss \jT « 
10 nm where T is the cavity finesse and AT w 50 mK. 
These values give the conditions on the length and tem- 
perature control loops to remain within the locking zone. 
These constraints arc compatible with the current perfor- 
mances of length and temperature controls. 

4 Linear cavity type II OPO 

In this section, we study the linear cavity case which is 
actually used in most experiments. We show here that the 
linear and ring cavity OPO have different behaviors when 
one takes into account the reflection phaseshifts on the 
cavity mirrors for the different interacting waves. 

One mirror, Mi is highly reflective for signal and idler 
and serves as a coupling mirror for the pump while the 
other mirror, M2 is highly reflective for the pump and 
serves as a coupling mirror for signal and idler (fig . 
The phase of the reflection coefficient for signal and idler 



are taken equal. 

We redefine the waveplate coupling constants since the 
signal and idler beams pass two times in the waveplate 1 : 

a = cos(A(f>) + i sin(Z\0) cos(2/}) = a e^ (10) 
e = isin(^)sin(2/3) (11) 

A<fr/2 being replaced by A<j>. 

As mentioned in the introduction and in the previous 
section, the linear cavity OPO has distinct features when 
compared to the ring cavity while the triple resonance does 
not change the behavior of the system. In the linear cavity 
the beams undergo two interactions per round-trip. As the 
phase is important in a parametric interaction the phase- 
shift between signal and idler and the pump beam between 
the two nonlinear interactions in the crystal must be taken 
into account. The equations for the field enveloppes at face 
(a) of the crystal can be written to the first order in g' : 



A 1 = a a r'e i{ -*-^[A 1 



u r (1 + e*)g'A At] 
[A 2 + (1 - e^g'AoA*,} 
A 2 = aor'eW) [A 2 + (1 + e*)cfMA\] 



+er ' e iS\ 



(12) 



+er'e l5 [A 1 + {\-e^)gA A* 2 



where 



UJ Q 



S = -^{2ne + 2nl + 2L) + Ci + C2 
2c 

S' = 0-i/) 

CJ 



<5o 



— (2n / + 2L) + Co 
c 

UJQ 



— (27i ~ 2n)l ~ —n(2e) + Co - 2C 2 
2c c 



(13) 
(14) 
(15) 
(16) 



2L is the total round-trip free propagation length, n and 
9 have been defined in section |3J 

One sees on the first two equations of expression 1121 
that when the phase-shift £ is taken equal to 0, the equa- 
tions are similar to the ring cavity case 2 , but with a crys- 
tal of double length (factor 2g'). This is no longer the case 



the free propagation and reflection on the coupling mirror 

simply shift the two waves by the same phase which does not 

change the effect of the waveplate 
2 when one neglects the second order term in eg if p is taken 

to be small 
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Fig. 7. Locking zone as a function of cavity length and crystal 
temperature for two values of £ : top £ = 0, bottom £ = 7r/4. 
The other values are the same : p = 5°, a = 3. 

when this parameter is changed. A non-zero value of £ has 
been shown to increase the threshold of a standard OPO 
by a significant amount [7j. In the case of a linear cavity 
with a birefringent element, a dissymmetry appears in the 
equations due to the terms 1 ± e^. Figd shows an exam- 
ple of the results obtained : one notices the dissymmetry 
between the two locking zones. 

Fig.[S]presents the value of the threshold on resonance, 
a res as a function of the temperature ST and the phase- 
shift £ for p — 5°. One observes that this threshold is 
no longer obtained for ST = as is the case for a ring 
cavity. For small values of p and £, a res remains reason- 
able (cr res < 3) inside a temperature range of approxi- 
mately 1 K. This value is small compared to the pure 
phase matching temperature range of 15 K. However, as 
£ increases and goes to n, a res diverges : as £ is fixed by the 
exact mirror dielectric structure, it is not adjustable ex- 
perimentally (except by changing the mirrors) : this can 
be a severe limitation to operation of the phase-locked 
OPO for small values of the waveplate angle p. 

When p is increased, the locking zone size increases 
and a larger range of temperature can be used with a low 
threshold. Figure shows the behavior of the normalized 
threshold on resonance a res as a function of £ and ST 
for p = 45°. In this case, the minimum value of a res is 
obtained for £ = it and ST — 0. When £ is lowered to 0, 
a res increases. The maximum value of a res is obtained for 
£ = and two values of the temperature : it is equal to 



. (part 1) : classical operation 7 
ST(K) 




Fig. 8. Normalized threshold on resonance oy es as a function 
of the crystal temperature 8T and of the phase-shift £ for p — 
5°. Unshaded surfaces correspond to values of (ST,SL) where 
frequency degenerate operation is not possible. 




Fig. 9. Normalized threshold on resonance o res as a function of 
the crystal temperature 8T and of the phase-shift £ for p — 45° 
(top). Same curve optimized for 8T (bottom) 

1.92 times the standard OPO threshold. This increase by 
a factor 1.92 is also found for £ = tt in the case of the 
standard OPO 0. 



8 



L. Longchambon et al: Non- linear and quantum optics of a type II OPO . . . (part 1) : classical operation 



5 Conclusion 

We have studied a system composed of an Optical Para- 
metric Oscillator containing a birefringent waveplate in- 
side the optical cavity. As shown previously 02) j this sys- 
tem allows phase locking of the signal and idler fields. 
We have obtained equations that are valid for all wave 
plate angles as well as in different cavity configurations, 
namely ring or linear cavities. We have shown that the 
zone where phase locking occurs can be described by the 
cavity length and the crystal temperature and consists of 
two zones. As the waveplate angle is increased, the size of 
the locking zone increases. The optimal configuration is 
obtained by inserting a A/2 waveplate in a ring cavity or 
a A/4 waveplate in a linear cavity with a 45° angle with 
respect to the crystal's axis. In the case of a ring cav- 
ity, the minimum threshold is obtained for a temperature 
such that the crystal birefringence compensates all the 
other birefringence in the cavity (waveplate and mirrors) 
and is equal to the standard OPO threshold. The effect 
of phase mismatch between the three waves is small for 
small values of the waveplate angle since the locking zone 
extension in temperature is small. As p is increased, the 
effect of phase mismatch becomes noticeable and limits in 
practice the extension of the locking zone. In a linear cav- 
ity, the mirrors phase-shift modifies the minimum thresh- 
old which becomes dependent on the waveplate angle and 
can become twice as large as the standard OPO thresh- 
old. This increase is known even in standard OPOs but 
a linear cavity is usually chosen for experimental reasons 
(losses, mechanical stability. . . ). In both cases (standard 
and self-phase-locked OPO), this increase is accompanied 
of a shift in the optimal crystal temperature which may 
be large and must be taken into account to operate the 
OPO at low threshold. 
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Appendix 

We give here the exact expression for the lower oscillation 
threshold in the case of a ring cavity : 

u — \/v 



I 



in 



(17) 



with 

u = ( 2 + r' 2 — 2r'ao cos(S) cos 
+a 2 cos (6 - 2ip) 



2vb 



(18) 



r' 2 + £o ~ 2r ' a o cos(<5) cos ( ^ - ip ) + a 2 , cos 2 (6> - 2tp) 



-1 - r' 4 - 2r' 2 al - 2r' I r' cos(2<5) + a x 



2tp) -2(1 + r' 2 ) cos((5) cos ( - - i/) 



r ao cos( 

with the parameters defined in the text. 
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